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ABSTRACT 

The dense concentration of stars and high velocity dispersions in the Galactic centre 
imply that stellar collisions frequently occur. Stellar collisions could therefore result 
in significant mass loss rates. We calculate the amount of stellar mass lost due to 
indirect and direct stellar collisions and find its dependence on the present-day mass 
function of stars. We find that the total mass loss rate in the Galactic centre due to 
stellar collisions is sensitive to the present-day mass function adopted. We use the 
observed x-ray luminosity in the Galactic centre to preclude any present-day mass 
functions that result in mass loss rates > 10~^MQyr~^ in the vicinity of ^ 1". For 
present-day mass functions of the form, dN/dM oc M~", we constrain the present-day 
mass function to have a minimum stellar mass < TMq and a povifer law slope > 1.25. 
We also use this result to constrain the initial mass function in the Galactic centre by 
considering different star formation scenarios. 
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1 INTRODUCTION 



The dense stellar core at the Galactic centre h as a radius of 
^ 0.15-0.4pc, a stellar d e nsity > lO'^Mppc" ^ llGenzel et al 
19961 : lEckart et all 1 19931 : iGenzel et aLll2003l : ISchodel et al 
2003) high velocity dispersions 100km s"'^), and Sgr A* 



the central supermassive black hole with a mass « 4 x IO'^Mq 
llEckart et all I2OO2I : ISchodel et al] |2002| . l2003l : iGhez et all 
12001 1 20081 ). Due to the extreme number densities and ve- 
locities, stellar collisions are believed to play an important 
role in shaping the stellar structure around the Galactic cen- 
tre, an d in disrupting the evolution of its stars. iGenzel et al.l 
l|l996l) found a paucity of the brightest giants in the galac- 
tic center, and proposed that collisions with main sequence 
(MS) stars could b e the culpri t. This hypothesis was found 
to be plausible by lAlexandeil (flgga ). Other investigations 
of collision s between giants and MS, white dwarf and neu- 
tron stars (|Bailev fc Daviej|l999l ) and c ollisions between gi - 
ants and binary MS and neutron stars (jPavies et al.|[l998l ) 
could not account for the dearth of o bserved giants. Th e con- 
tradictory results were resolved bv lDale et al] l|2009i ). who 
concluded that the lack of the faintest giants (but not the 
brightest giants) could be explained by collisions between gi- 
ants and stellar mass black holes. Significant mass loss in the 
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giants' envelopes after a collision would prevent the giants 
from becoming bright enough to be observed. 

The above studies concentrated on collisions involv- 
ing particular stellar species with particular stellar masses. 
To examine the cumulative effect of collisions amongst an 
entire ensemble of a stellar species with a spectrum of 
masses, one must specify the present-day stellar mass func- 
tion (PDMF) for that species. The PDMF gives the current 
number of stars per unit stellar mass up to a normaliza- 
tion constant. Given a certain star formation history, the 
PDMF can be used to determine the initial mass function 
of stars (IMF) , the mass function with which the stars were 
born. There is currently no consensus as to whether the 
IMF in the Galactic centre devi ates from the canonical IMF 
jBastian. Covev. fc Meve3l201Gl ). 

First described by Salpeter more than 50 years ago 
l|Salpeteij[l955l ). the canonical IMF is an empirical func - 
tion which has been found to be universal (|Kroupall200ll ). 
with the Galactic centre as perhaps the sole exception. 
iManess et alJ (|2007l ) found that models with a top-heavy 
IMF were most consisten t with observat ions of the cen- 
tral parse c of the Galaxv. [Paumard et al.l (|2006l ). and sub- 
sequently [fiartko^^eF^l] (|2010l ) found observational evidence 
for a flat IMF for the y oung OB-stars in the Galactic cen- 
tre. O n the other hand, iLockmann. Baumgardt. fc Kroupal 
l|2010l ) concluded that models of constant star formation 
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with a canonical IMF could explain observations of the 
Galactic centre. 

In this work we use calculated mass loss rates due to 
stellar collisions as a method to constrain the PDMF for 
main sequence stars in the Galactic centre. We construct a 
simple model to estimate the actual mass loss rate in the 
Galactic centre based on observed x-ray emission. PDMFs 
that predict mass loss rates from stellar collisions greater 
than the observed rate are precluded. This method allows 
us to place conservative constraints on the PDMF, because 
we do not include the contribution to the ma ss loss rate from 
stella r winds from massive evolved stars (iBaganoff et al.l 
I2OO3I ). Specifically, this method allows us to place a lower 
limit on the power-law slope and an upper limit on the 
minimum stellar mass of the PDMF in the Galactic centre 
(see § [SJ . Inclusion of the mass loss rate from stellar winds 
(or other sources) could further constrain the PDMF of the 
Galactic centre. 

We present novel, analytical models to calculate the 
amount of stellar mass lost due to stellar collisions between 
main sequence stars in § [2] through § 12.31 In § [3] we develop 
the formalism for calculating collision rates in the Galactic 
center. We utilize our calculations of the mass loss per col- 
lision, and the collision rate as a function of Galactic radius 
to find the radial profile of the mass loss rate in § |4l Since 
the amount of mass lost is dependent on the masses of the 
colliding stars, the mass loss rate in the Galactic centre is 
sensitive to the underlying PDMF. By comparing our calcu- 
lations to mass loss rates obtained from the x-ray luminosity 
measured by Chandra, in §[5] we constrain the PDMF of the 
Galactic centre. We derive analytic solutions of the PDMF 
as a function of an adopted IMF for different star formation 
scenarios, which allows us to place constraints on the IMF 
in §6. In § 7, we estimate the contribution to the mass loss 
rate from collisions involving red giant (RG) stars. 



2 CONDITION FOR MASS LOSS 

Throughout this paper we refer to the star that loses mate- 
rial as the perturbed star, and the star that causes material 
to be lost as the perturber star. Quantities with the subscript 
or superscript "pd" or "pr" refer to the perturbed star and 
perturber star respectively. We work in units where mass is 
measured in the mass of the perturbed star, Mpd, distance in 
the radius of the perturbed star, rpd, velocity in the escape 
velocity of the perturbed star, v^^f^ (= \/'2GMpd/rpd), and 
time in rpd/v^f^. We denote normalization by these quan- 
tities (or the appropriate combination of these quantities) 



^ Note that for any particular collision, it is arbitrary which star 
we consider the perturber star, and which star the perturbed star. 
Both stars will lose mass due to the presence of the other, so in 
order to calculate the total mass loss, we interchange the labels 
(pd-f->pr), and repeat the calculation. 
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(1) 

We refer to collisions in which b > Vpd + Vpr as "indirect" 
collisions, and collisions in which b ^ rpd + rpr as "direct" 
collisions. The impact parameter, b, is the distance of closest 
approach measured from the centers of both stars. 

We consider the condition for mass loss at a position, 
f, within the perturbed star to be that the kick velocity 
due to the encounter at f exceeds the escape velocity of the 
perturber star at f, A'3(f) ^ ■Desc(r). The escape velocity 
as a function of position within the perturbed star can be 
found from the initial kinetic and potential energies of a test 
particle at position f , 



+ 47r / p(r')r'dr' 



(2) 



where Mint is the mass interior at position f ' and p is the 
density profile of the star. 



2.1 Mass loss due to Indirect Collisions 

To calculate the mass lost due to an indirect collision, we 
first calculate the kick velocity given to the perturbed star 
as a function of position wi thin the star . We work under 
the impulse approximation (|Spitzerlll958l ). valid under the 
condition that the encounter time is much shorter than the 
characteristic crossing time of a constituent of the perturbed 
system. 

Given a mass distribution for the perturbed system, ppd 
and a potential for the perturber system, $, the kick velocity 
aft er an encounter under the i mpulse approximation is given 
bv lBinnev fc Tremaind lj2008l 'l: 



Aw(r) 



1 



Ppd(P,t)V$(P,t)dV] dt. 



(3) 



Equation @ can be simplified by expanding the gradient of 
the potential in a Taylor series, resulting in 



Au (r) 



2GM, 




+ 



(4) 



The expansion is valid under the "distant tide" approxima- 
tion which is satisfied when rpd <^ b. The parameter Vrei is 
the relative speed between the stars {vrei = \vpd — «pr|). We 
are interested in the magnitude of equation Q , which when 
normalized to the units that we have adopted for this paper 
is 



Av{x, y) = 7-\/ £2 + 



(5) 
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Table 1. Coefficients of polynomial fits for A (7) with varying 
polytropic indeces 



n 


at) 


ai 


02 


13 


1.5 


0.395 


-0.865 


0.559 


-0.091 


2.0 


0.210 


-0.424 


0.246 


-0.032 


2.5 


0.105 


-0.197 


0.102 


-0.101 


3.0 


0.051 


-0.088 


0.040 


-0.003 



where 



(6) 



To solve for the mass lost per encounter as a function 
of 7, we consider a star within a cubic array, where the 
star contains ~ 3 x 10® cubic elements. As a function of 7 
we compare the kick velocity in each element to the escape 
velocity for that element, and consider the mass within the 
element to be lost to the star if the velocities satisfy the 
condition given § (2] We note that by ^ 10^ elements, the 
results converge to within about 2%, and we are therefore 
confident that ~ 3 x 10® provides adequate resolution. 

To calculate the amount of mass in each element, the 
density profile for the perturbed star must be specified. As 
with s everal previous studies on mass loss due to stellar coUi- 
sions teenz fc Hillsl l 19871 . [lOQ^ : iLai. Rasio. fc Shapiroll 19931 : 
[Rauch 1999) wo utilize polytropic stellar profiles. Polytropic 
profiles are easy to calculate, and yield reliable results for 
stars of certain masses. Polytropic profiles of polytropic in- 
dex n — 1.5 describe the density structure of fully convec- 
tive stars, and therefore very well describe MS stars with 
Mi, < 0.3Mq (nearly fully convective) and MS stars with 
M* > IOM0 (convective cores). MS stars with M, > IMq 
have radiative envelopes, and are therefore well described 
by n = 3. For n for stars with masses of 0.3 — IM0 and 
5 — IOMq, we linearly interpolate between n = 1.5 and 3. 
We discuss the uncertainties introduced by this approach in 
§111 Note that this approach is biased towards zero-age main 
sequence stars, since as stars evolve, they are less adequately 
described by polytropic profiles. 

We plot the fraction of mass lost from the perturbed 
star per event. A, as a function of 7 in Fig. [T] for several 
polytropic indeces. The lines are third order polynomial fits 
to our results, in the range of 0.98 ^7^5. We list the co- 
efficients of the polynomial fits in Table [T] For each density 
profile, no mass is lost up until 7 of about 0.98, and there- 
after the mass loss increases monotonically. The increasing 
trend is due to the fact that larger perturber masses and 
smaller impact parameters result in an increased potential 
felt by the perturbed star. Smaller velocities also cause more 
mass to be lost, as this increases the "interaction time" be- 
tween the perturber and perturbed stars. 

The location of the mass loss within the perturbed star 
for fixed 7 depends upon the polytropic index, since the 
escape velocity within the star is dependent upon the density 
profile, as indicated by equation. ([2]). In Fig. (2] we illustrate 
where mass will be lost in the perturbed star by plotting 
contours of the kick velocity {Av{f)) due to the encounter 
normalized to v^tdf) for n = 1.5 and n = 3 (top and bottom 
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Figure 1. The fraction of mass lost per collision as a function 
of 7 for several poljrtropic indeces. The lines are third order poly- 
nomial fits, whose coefficients are given in Table [Tl 



row respectively). We show two different cases: a slightly 
perturbing encounter with 7 = 1.2 in the first column, and 
a severely perturbing encounter with 7 = 1.6 in the second 
column. The grey region underneath shows where mass is 
still left after the encounter, since Av/v^f^{f) within this 
region is < 1. The 7 = 1.6 encounter results in bigger kick 
velocities, and so we see that the mass loss penetrates farther 
into the star. We note that the shape and magnitude of the 
contours for both polytropic indeces at fixed 7 converge at 
large radii. This is due to the fact that regardless of the 
polytropic index used, v^f^ converges to the same value at 
large radii when the second term in equation ([2]) becomes 
negligible. Even though the location of where mass is lost 
is similar for different polytropic stars at the same value of 
7, the amount of the mass lost is substantially different (as 
shown in Fig. [T]), due to the different density profiles. 

2.2 Validity of approach for indirect collisions 

The impulse approximation is valid provided that the time 
over which the encounter takes place, tenc, is much shorter 
than the time it takes to cross the perturbed system, t cross ■ 
To estimate when our calculations break down, we approxi- 
mate tenc as b/vrei, and tcross as ts, the time it takes for a 
sound wave to cross an object that is in hydrostatic equilib- 
rium: 

1 1 



ts 



These approximations lead to the condition that 



(7) 



t'rei" ^ (8) 

lAguilar fc Whitd (|l985l ) find that for a large range of 
collisions, the impulse approximations remains remarkably 
valid, even when t^nc is almost as long as tcross- We there- 
fore assume that the impulse approximation holds until the 
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Figure 2. Slices through the perturbed star along the plane parallel to the perturber star's trajectory. The first column (panels a) and 
c)) correspond to encounters with 7 = 1.2, and the second column (panels b) and d)) correspond to 7 = 1.6. The first row (panels a) 
and b)) have n = 3, and the second row (panels c) and d)) have n = 1.5. The contours are the kick velocity within the star due to the 
encounter normalized to the escape velocity (as a function of r). The outline of the grey region underneath has Av/v^fc{f) = 1, so that 
the grey region represents the location of where mass is still left after the event. 



left hand side of equation ([§1 is ~ 1. Our calculation of A 
as a function of 7 should therefore be valid for 7 < jvaUd, 
where jvaiid = Mpr/b^- We plot contours of log(7„aHd) in the 
Mpr/Mpd - b/rpd parameter space in Fig. [31 where both the 
X and y axes span ranges relevant to our calculations. The 
shaded grey area in the figure is the region of the parameter 
space where the impulse approximation predicts non-zero 
mass loss due to the encounter. The figure shows that "/vaud 
is smaller for low Alpr to Mpd ratios at high impact param- 
eters. In fact, most of the right side of the parameter space 
has "/valid less than 0.98 (where below this value, the impulse 
approximation predicts no mass lost). 

In our calculations, when, for any particular set of 
Mpr/Mpd and b/r-pd, 7 > jvaiid, we adopt A(7 > jvaiid) = 
A(7 = "/valid)- This approach represents a lower limit on 
the amount of mass loss that we calculate, since mass loss 
should increase with increasing 7. We find, however, that if 
we set A(7 > "/valid) = 1 (which represents the absolute up- 
per limit in the amount of mass lost) the change in our final 
results is negligible at small Galactic radii. At large radii. 



where the mass loss from indirect collisions dominates (see 
§ 13} , the results change by at most a fact of ~ 2. 



Equation ^ was derived under the assumption that 
the impact parameter is much bigger than both rpd and rpr- 
Since Aw scales as b~^ , the equation predicts that most mass 
loss occurs for small impact parameters. However, given the 
assumption that was used to derive the equation, the regime 
of small impact parameters is precisely where equati on. Q 
breaks down. Numerical simulations jA guilar fc White! 
'l985!; iGnedin. Hernguist. fc Ostrikeij [19991 ) show that for a 
variety of perturber mass distributions, the energy input into 
the perturbed system is well described by equation Q for 
b ^ Sr^, where is the half mass radius of the perturber 
system. For an n = 3 polytropic star, br^ = 1.4r*. Since for 
indirect collisions, b/rpd = 1 + rpr/rpd + d/rpd (where d is the 
distance between the surface of both stars), there is only a 
small region in our calculations, ^ d/rpd < (OA — rpr/rpd), 
for which equation may give unreliable results. 
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Figure 3. Contours of log{-y„alid) in the Mpr/Mp^i - hjr^^ pa- 
rameter space, where 'y^alid is defined in § [2] The shaded grey 
area indicates where the impulse approximation predicts non-zero 
mass loss 



log(M/MJ 

Figure 4. Mass-radius relations used in studies of calculating 
mass loss from stars due to stellar collisions. The thin lines are 
power-law relatio ns of power-law index 1.0, .8, a nd 0.85 used 
by, Ranch (;1999| ). iLai. Rasio. fc Shapird l|l993l ) and lBenz fc Hillj 
lll992l) respecti v ely. T he dotted line is the relation used by 
iFreitag fc Bend l|2005h . and the thick line is the relation used 
in this work. 



2.3 Mass loss due to direct Collisions 

A number of papers over the past few decades have ad- 
dressed the outcomes of stellar collisions where the two stars 
come so close to each other that not only gravitational, 
but also hydrodynamic forces must be accounted for. Early 
studies used one or two dimensional low resolution hydrody- 
namic simulations (e.g. Mathis 1967; DeYoung 1968). Mod- 
ern studies typically utilize smooth particle hydrodynamics 
with various stellar models, mass-r adius relations, and vary- 
ing degrees of particle resolution (iBenz fc HilldligSTl . 1 19921 : 

■ i 



iLai. Rasio. fc Shapirc|[l993l : [ Rauc hl l 19991 ). A detailed review 



of the literature can b e found in this area can be found in 
iFreitag fc Ben3 l|2005h . 

We approach the problem of direct collisions in a highly 
simplified, analytic manner without hydrodynamic consid- 
erations, and find that for determining the amount of mass 
lost, our method compares well to the complex hydrody- 
namic simulations. As a first order model, we approximate 
the encounter as two colliding disks, by projecting the mass 
of both stars on a plane perpendicular to the trajectory of 
the perturber star. The problem of calculating mass loss 
then becomes easier to handle, as it is two dimensional. We 
also assume that mass loss can only occur in the geometrical 
area of intersection of the two stars. 

We find the kick velocity as a function of position in 
the area of intersection by conserving momenta, and by as- 
suming that all of the momentum in the perturber star in 
each area element was transfered to the corresponding area 
element in the perturbed star. Working in the frame of the 
perturbed star, and with a polar coordinate system at its 
center (so that r = ^/a;^ + y^), we find 



Av{f) = 



Epd(r) 



(9) 



The parameters Sp^ and Tipd represent the surface density of 
the perturber and perturbed stars respectively (E = J pdz). 

To find the region of intersection, we need to know the 
impact parameter, and the radii of both stars. To obtain the 
stellar radii as a func tion of mass, we use the mass-radius 
relation calculated by iKippenhahn fc Weigerj ([l994) for a 
MS star with Z = Zq, Xh = 0.685, and X^e = 0.294 from 
a stellar evolution model, where X represents the mass frac- 
tion. We fit a polynomial to their Fig. 22.2, and extrapolate 
on the high and low mass ends so that we have a mass- 
radius relation that spans from about 0.01 to 150Mq. We 
compare our mass-radius relation to those used in other 
studies of direc stella r collis ions i n Fig. HI iRauch (|199^ , 
iLai. Rasio. fc Shapird (|l993l ) and iBenz fc Hilj (|l992l ) all 
adopted power laws wit h power law ind i ces of 1.0, 0.8, 0.85 
respectively (thin lines). IFreitag fc Bend l|2005l ) (dotted line) 
use main sequence stellar evolution codes to obtain a mass- 
radius relation for masses > 0.4Mq, and a polytropic mass- 
radius relation of n = 1.5 for masses < 0.4Mq. 

Our simple model for calculating mass loss due to di- 
rect stellar collisions compares surprisingly well to full blown 
smooth particle hydrodynamic simulations. We borrow plots 
of the fractional amount of mass lost as a function of impact 
para meter for specific relat ive velocities and stellar masses 
from lFreitag fc Bend l|2005l ) (Figs.[S]and[n}. They show their 
own work, the best calculations of mass loss due to stellar 
collisions to date. For comparison, and to show how the cal- 
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culations have evolved over the years, the results from older 
studies are also shown. Our own results are plotted (dashed- 
dotted black lines) over these previous studies. We make sure 
to show results spanning a wide range of stellar masses and 
relative velocities. Note that these plots show the fractional 
amount of mass lost from both stars normalized to the ini- 
tial meisses of both stars, and that the impact parameter 
is normalized to the sum of both stellar r adii. Our results 
show t he same qualitative trends seen in the lFreitag fc Bend 
curves, even replicating several "bumps" seen in their 
cur ves (see the two bottom panels of Fig.|6]). As compared to 
the lFreitag fc Bend (|2005l ) results, for any specific set stellar 
masses, relative velocity and impact parameter, our calcula- 
tions sometimes over or under-predict the amount of mass 
lost by of a factor of a few to at most a factor of 10. We 
discuss the error introduced into our main calculations by 
this discrepancy at the end of ^ 



3 STELLAR COLLISION RATES IN THE 
GALACTIC CENTER 

To calculate mass loss rates in the Galactic center, we will 
need to find the collision rates as a function of the perturber 
and perturbed star masses, impact parameter, and relative 
velocity. Additionally, the collision rate will be a function of 
distance from the Galactic center, since the stellar densities 
and relative velocities vary with this distance. In this section, 
we first present the Galactic density profile that we use, and 
we then derive the differential collision rate as a function of 

these parameters. 

W e adopt the stellar density profile of ISchodel et al.l 
one of the best measurements of the density pro- 
file within the Galactic centre to date. Using stellar counts 
from high resolution images of the galactic center, they find 
that the density profile is well approximated by a broken 
power law. Moreover, they use measured velocity disper- 
sions to constrain the amount of enclosed stellar mass as a 
function of galactic radius, rgai- Using their density profile, 
and velocity dispersion measurements, they find that 

2.8 ± 1.3 X 106Mopc-3 (^J^y^'^ for r^ai < 0.22 pc 
2.8 ± 1.3 X 10*^ (^J^y^'^'' for r^ai > 0.22 pc 

(10) 

Their average density can be converted into a local density, 
piXgai), by considering the definition of p, 

4/3^^3^^ > (11) 

from which we derive: 

p(r,.0 = P(^«.0 + ^^9^. (12) 

We use equations (|10|) and (|12|l . to find p{r), and plot 
the result in Fig. [T] We "smoothed" the unphysical discon- 
tinuity in p arising from the kink of the broken power law 
fit by fitting a polynomial to equation (|10|l . 

The differential collision rate, dV, between two species. 




Figure 7. The stellar dens ity profile that we a dopt, based on 
the average density profile of lSchodel et al.] l l2007t) . 

"1" and "2" at impact parameter b characterized by distri- 
bution functions /i and /2 , and moving with relative velocity 
l«i — i^bl in a spherically symmetric system is 

dF = fl{rgaUVl) d'^Vif2{rgal,V2) d^V2 

X \vi - V2\2nbdb4TTrgaidrgai- (13) 
For simplicity, we adopt Maxwellian distributions, 

r r ^ X ni.2{rgal) -„2 /2ct2 

fi,2{r,ai,vi,2) = (27ra2)3/2 ^ ■ ' (14) 

where we find the velocity dispersion, a, from the Jean's 
equations. Assuming an isotropic velocity dispersion, a 
spherical distribution of stars, and a power-law density pro- 
file with power-law slope 13, p oc r~fi, the Jean's equations 
lead to cr^ = r^f^+F) ' Msmbh = 4 x IO^Mq. From 

Fig. [71 it is evident that P{rgai), but for simplicity, we adopt 
an averaged value of /3, /3 = 1.3. Note that we have also as- 
sumed that the enclosed mass at position Vgat is dominated 
by the SMBH. This assumptions is valid out till ~ Ipc, which 
is also the point where our impulse approximation starts to 
break down. 

A change of variables allows one to integrate out 3 of 
the velocity dimensions, and to write the expression in terms 
of Vrei (see Binney & Tremaine 2008). We can also take 
into account the fact that both species have a distribution 
of masses by introducing, ^1,2, the PDMF, which gives the 
number density of stars per mass bin (^ = dn/dM). We 
adopt a power law PDMF, 

^ocAf-", (15) 

that runs from some minimum mass, Mmin to a maximum 
mass Mmax- Since most initial mass functions are param- 
eterized with a power law, the present-day mass function 
might be modified from a power law due to the effects of 
collisions and stellar evolution. Although the actual PDMF 





Figure 5. The calculated fractional amount of mass lost as a function of impact paramet er from several works. Our res u lts ar e 
the black dashed-dotted lines. The acroiiyms F B05, R99, LRS93, BH87, a nd SS66 refer to iFreitag fc Bend ll2005h . iRauchI lll999^, 
iLai. Rasio. fc Shapirol lll993l ). rBen2 &: Hillsl lll987fl and lSpitzer fc Saslawl lll966l ) respectively. The figures are adopted from lFreitag &: Bend 
(|2005l '). 



might have deviations from a power law, adopting a power 
law provides us with a quick and simple way to parameter- 
ize the PDMF. Taking all of this into account, and assuming 
that the relative velocities are isotropic, we arrive at the final 
non-dimensionalized expression for the differential collision 



rate: 

X Mi°'M2"flalbdbdfgal dVrel dMi dMi- (16) 

The tildes denote normalization by the proper combination 




M.=0.5 Mq M„=2.0 Mq V~/V^ = 0.Q7 



0.1 



0.01 



10- 



10- 




0.2 0.4 
b/{R,+R„ 



0.1 



0.01 



+ 



10-3 ^ 



10-'' ^ 



10- 



T 




0.2 



0.4 
b/{R,+R,) 



0.6 



Figure 6. The calculated fractional amount of mass lost as a function of im pact parameter from several works. Our results are the 
black dashcd-dotted lines. The figures are adopted from .Freitag fc Ben3 l|2005l ). 



of M2, r2, and Uesc- The parameter K is the normalization profile of Fig. [7| and the following expression: 
constant for ^, which can be solved for by using the density rMmam j„ 

^MdM 
dM 



^j^(raaO o • (17) 



Constraining the Stellar Mass Function in the Galactic Center Via Mass Loss from Stellar Collisions 



Since the expression for K, whicii controls the total num- 
ber of stars, has no time dependence, our expression for the 
PDMF assumes a constant star formation rate in the Galac- 
tic center. 

Our calculations involve the computation of multidi- 
mensional integrals over a two dimensional parameter space 
(see 311). Therefore, for the ease of calculation, we ignore the 
enhancement of the collision rate due to the effects of grav- 
itational focusing. This results in a conservative estimate of 
the collision rate. As two projectiles collide with each other, 
their mutual gravitational attraction pulls them together, 
resulting in an enhancement of the cross section: 



S{1 + 



2G(hh + M2 

ret 



(18) 



We discuss the uncertainties in our final results due to ig- 
noring gravitational focusing at the end of 3S1 

To illustrate the frequency of collisions in the Galactic 
center, we integrate equation. (|16p over Vre^i, Mi and M2 as- 
suming a Salpeter-like mass function (a = 2.35, Mmin = 
O.IM0 and M^ax = 125M©) to obtain dV /{dlnrgaidb) 
as a function of rgai (Fig. [H^fl We plot dV /{dlnrgaidb) 
for several different impact parameter values. We calculate 
dV / (dlnrgaidb) with and without the effect of gravitational 
focusing (solid and dashed lines respectively). The latter is 
obtained by multiplying equation (|16p by the gravitational 
focusing enhancement term before the integration. As ex- 
pected, gravitational focusing is negligible at small Galactic 
radii since typical stellar encounters involve high relative 
velocities. As the typical relative velocities decrease with 
increasing Galactic radius, the enhancement to the colli- 
sion rate from gravitational focusing becomes important. 
The figure also shows that gravitational focusing becomes 
less important with increasing impact parameter since the 
gravitational attraction between the stars is weaker. Fig. [HJj 
shows the cumulative differential collision rate (integrated 
over rgai) per impact parameter as a function of rgai- Again, 
we plot the results with and without gravitational focusing 
and for the same impact parameters. 



4 MASS LOSS RATES IN THE GALAXY 

To calculate the mass loss rate from stars due to collisions 
within the Galactic center, we multiply equation (jlGp by the 
fraction of mass lost per collision, A (7), and compute the 
multi-dimensional integral. We calculate the total mass loss 
rate from both the perturbed and perturber stars by simply 
interchanging the "pr" and "pd" labels and re-performing 
the calculation. 

We first compute the differential mass loss rate for in- 



^ Tills figure, and subsequent figures in tliis paper with rgai 
the independent variable start from rgai = 10~®pc. This value of 
rgai corresponds to the tidal radius for a IMq star associated with 
a 4 X 1O^M0 SMBH. Although stars of different masses will have 
slightly different tidal radii, the main conclusions of our paper arc 
based off of distances in rgai of order O.lpc (see [jSj, well above 
the tidal radius for any particular star. 



direct collisions. The mass loss per collision is given by: 

( for -y < 0.98 

A(7) = I polynomial for 0.98 < 7 ^ ivaiid ■ (19) 

[ /^{■~1vaUd) for 7 > 7„aHd 

The coefficients for the polynomial depend on the polytopic 
index of the perturbed star (and thus on its mass) and are 
taken from Table [T] We multiply equation (|19|) and equa- 
tion H16p and simplify the integration. In principle, b should 
go to 00, but we cut off the integral at bmax ~ 20 as we 
find that the results converge well before this point. The ve- 
locity integral is also cut off at Vmax due to the fact that 
A(7) becomes zero below 7 = 0.98. This cut-off corresponds 
to Vmax ~ ^^^ggj^a'"'"' ■ We may safely throw away the expo- 
nential as Vmax ^ S-'^{fgai) for the range of Pgai that we 
consider. Thus, the integral that we evaluate is: 



dM 
dlnfgai 



^ix^l'^b-^flaiK' 



pd 



rbmax 



X / bv:^,iApdi'y)Mpr"Mpj;'dbdMpr 

X dMpadVrei- (20) 

For direct collisions, A{b, Mpr, Mpd,Vrei) is calculated 
given the prescription in Sec l2.3l To evaluate the multidi- 
mensional integral, we make the approximation of evaluat- 
ing Apd at Vrei = 2(T. The factor of Apd thus comes out of 
the Vrei integral, so that the Vrei integral can be performed 
analytically: 



dM 



dlnVgal I JM^i^ J M,^.^ Jo 

/ pa IILLIL IILlll, 

X bApd(b, Mpr, Mpd, Vrei = 2a{fgal)) 

X Mpr^M-fdbdMprdMpd. (21) 

We evaluate the remaining integrals numerically. 

Once values for a, Mmin and Mmax are specified, equa- 
tions (|20|l and (I21|l can be integrated to obtain the mass 
loss rate as a function of Galactic radius. To show how the 
mass loss rate profiles vary with Mmin, Mmax and a, we plot 
dM / dlnr gal for direct collisions in Fig. [5]and vary these pa- 
rameters. In the figure, we have evaluated Mmin at 0.05, 0.5 
and 5Mq, Mmax at 75, 100, 125Mq and a from 1.00 to 2.5 in 
equal increments. The parameter Mmin increases vertically 
from the bottom panel to the top, Mmax increases horizon- 
tally from the left panel to the right, and in each panel a 
increases from the bottom to the top. We have indicated a 
Salpeter-like mass function (a = 2.29, Mmin ~ 0.5Mq and 
Mmax = 125Mq) with the dashed line. Mass loss is extensive 
and approximately constant until about rgai of 10~'^pc and 
then drops dramatically. This drop reflects that fact that 
collisions are less frequent at larger radii since star densi- 
ties and relative velocities drop. The amount of mass lost 
for any direct collision also decreases with galactic radius 
since A decreases with decreasing relative velocities. Note 
that the proflles are approximately constant as a function of 



M„ 



so that the choice of Mmin determines the extent of 



the mass loss rate. 

In Fig. [To] we show the contributions to dM/dlnrgai 
from both direct and indirect collisions for A/min=O.2M0, 
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Figure 8. a) The differential collision rate per logarithmic Galactic radius per impact parameter as a function of Galactic radius for 
several different impact parameters. The solid (dashed) lines were calculated ignoring (including) gravitational focusing. The curves 
were made by made by integrating equation II16I I (with and without the gravitational focusing term) assuming Salpeter values, b) The 
cumulative differential collision rate (integrated over Vg^i) per impact parameter with and without gravitational focusing for the same 
impact parameter values, a and b) The vertical line in each panel is placed at rg„; = 0.06pc, the upper bound in our integration across 
^gal 3S performed in 



Mmax=WOM.Q and a — 1.2. We find that at small radii the 
mass loss rate is dominated by direct collisions, and at large 
radii it is dominated by indirect collisions. Mass loss due to 
indirect collisions is suppressed in the Galactic center, due to 
the very fast relative stellar velocities. Even though the high 
velocities (and high densities) in the Galactic center make 
collisions more frequent, under the impulse approximation, 
when velocities are very fast, mass loss is minimized. 

To illustrate which mass stars contribute the most to 
the total mass loss rate, we plot dM / dlnMpd as a function 
of Mpd in Fig. [TT] for several different PDMFs. The range 
of integration we choose for Vgai is from to 0.06pc (see 

We choose Mmin to be 0.05, 0.5 and 5MQ(left to right 
in the figure), and we use a constant Mmax of 125M0. In 
each panel, we vary a from 1.5 to 2.5 in equal increments. 
The figure shows that for Mmin = O.O5M0, changing a has 
little effect on what mass stars contribute the most to the 
mass loss rate (although, the total mass loss rate is decreased 
with increasing a). For the Mmin ~ O.5M0 and 5M0 cases, 
increasing a results in lower mass stars contributing more 
to the mass loss rate. This trend makes sense, since PDMF 
profiles with higher values of a have fractionally more lower 
mass stars. 

To test how our interpolation between the n = 1.5 and 
3 polytropic indeces affects the main results of this paper, 
we consider two extreme cases. The first case we consider 
has n = 1.5 for M^, < IMq and M^, > 5M0, and n = 3 for 
IM0 ^ M* ^ 5M0. This approach has n — 1.5 for much 
of the mass spectrum, and should result in the highest mass 
loss rates since (as is evident from Fig. [l| collisions with 
the perturbed star having n = 1.5 result in the most mass 
lost. This is due to the fact that for n — 1.5 stars, the mass 
is less centrally concentrated, and more mass can therefore 
escape at large radii which receive a stronger velocity kick. 



The second case we consider has n = 1.5 for < O.3M0 
and > IOM0, and n = 3 for O.3M0 ^ ^ IOM0. 
This case should result in the smallest mass loss rates, since 
it has n = 1.5 for a smaller fraction of the mass spectrum. 
Since different mass functions have different fractions of the 
total mass in the neighborhood of IM0 (where we expect 
the least mass loss per collision since n = 3), we test the 
two cases for several different mass functions. We find that 
differences in dM / dlnrgai{rgai) for both cases are relatively 
minor, differing at most by ~ 10% depending on the mass 
function that we use. 



5 CONSTRAINING THE MASS FUNCTION IN 
THE GALACTIC CENTER 

It is known through x-ray observations from Chandra, that 
the central supermassive black hole in the Galactic center is 
surrounded by gas donated from stellar winds (e.g. Baganoff 
2003). The x-ray luminosity is due to Bremsstrahlung 
emission from unbound material supplied at a rate of ~ 
lO"^M0yr"^ (|Naiarro et al.lll997l '). This unbound ra aterial 
has been studied theoretically by lQuataer"3 (|2004l ). who 
solved the equations of hydrodynamics (under spherical 
symmetry) to f ollow how the gas is accreted onto Sgr A*. 
iQuataertl (|2004l ) finds that his model agrees with the level 
of diffuse x-ray emission measured by Chandra, and predicts 
an inflow of mass at rgai ~ 1" at a rate of ~ lO~^M0yr~"'^. 

It is interesting to note that the mass loss rate profiles 
due to stellar collisions are consistent with the observed x- 
ray surface brightnes s profile. F r om hi s derived temperature 
and density profiles, IQuataertl (|2004l ) solves for the x-ray 
surface brightness as a function of rj„; . IQuataertl (120041 ) has 
adopted a source term, q{rgai) (the stellar mass loss rate 
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Figure 9. Mass loss rates as a function of Galactic radius due to direct collisions for various parameters of Mmin Mmax, and a. The 
parameter Mmin varies in each panel from bottom to top, and Mmax varies from left to right. The power law slope, a varies within each 
panel from 1.00 (top line) to 2.5 (bottom line) in equal increments of 0.188. The dashed line corresponds to a Salpeter-like mass function 
values (Mmin = 0.5Mq, Mmax = 125M0 , a = 2.29). The arrows indicate the range in the x-ray observations (rg^i < 1.5") with which 
we use to constrain the PDMF (see ifs) . 
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Figure 11. The amount of mass loss per logarithmic mass interval of the perturbed star as a function of the perturbed star's mass. 
Each line was calculated with a different PDMF. The titles in each panel indicate the value of Mmin used for that panel. In each panel, 
a goes from 1.5 to 2.5 in even increments of 0.167, and for each line Mmax = 125M0. 
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Figure 10. Mass loss rates due to direct and indirect stellar 
collisions within the Galactic center for Mmin=0.2MQ, Mmax = 
IOOM0 and a = 1.2. The arrow indicates the range in the x-ray 
observations {rg^i < 1.5") with which we use to constrain the 
PDMF (see JS}. 



per unit volume) of q{rgai) oc r^Jj for Tgai G [2", 10"], for 
T] = Q, 2 and 3. He finds that the model with ?7 = 3 closely 
fits Chandra's observed surface brightness profile. From the 
definition of q{rgai), dM/dlnVgai oc rl^,q{rgai) oc r^^'P . It is 
eviden t from Fig. [9] that in the range of rgai that lQuataer"3 
l|2004h considers (about 0.08 to 0.39pc), dM/dlnrgai is a 
power law. Further, depending on which particular set of 
mass function parameters that we consider, rj = 2.6 — 2.8 
for stellar collisions. These values of rj are very close to the 
model most consistent with the Chandra surface brightness 
profile. 

Using the the 2-lOkeV luminosity as measured by Chan- 
dra l|Baganoff et al.1 120031 ) . we estimate the total mass loss 
rate at a radius of rgai ~ 1.5" (0.06pc). We use the word 
"total" to indicate the mass loss rate integrated over Galac- 
tic radius. By using this total mass loss rate as an upper 
limit, we will be able to constrain the PDMF in the Galac- 
tic center by precluding any PDMFs with total mass loss 
rates greater than this value. We will do this by integrating 
our calculated mass loss rate profiles (e.g.. Figs. l9l and [10)) 
over rgai ■ 

Unbound material at a radius rgai has a dynamical 
timescale of 



idyn (^gai ) 



Vchar ^r gal ) 



1.1 X loVrs ( ^ 
pc 



(22) 



where the characteristic velocity at radius rgai, Vchar(rgai), is 
taken as the velocity dispersion as given in §[3] The electron 
density at radius rgai may therefore be estimated by: 

r \ r \ Mtdynifgal) 

neirgai) ~ npirgai) j — 3 



= 1.1 X lO^cm"^ 



M 



Moyr-i 



^ (23) 
pc ' 



where rrip is the proton mass. 



For thermal Brems strahlung emission, the volu me emis- 
sivity {dE/dVdtdv) is l|Rvbicki fc Lightmanlll979l ) 



JJ 



1.8 X 10"^*ergs"^cm"^Hz"^ 



2 fx 
K 



-1/2 



X e ' ^ gff, 



(24) 



where we set gff — 1. The luminosity in the 2-lOkeV band, 
L2-10, can be found substituting equation (|23|l into equa- 
tion (|24p and integrating the volume emissivity over volume 
(assuming spherical symmetry) and frequency: 



6.7 X lO^^ergs^^ 



M 



Moyr- 



pc 



^/2 ' 



VkeVy 



(25) 



We have assumed a constant temperature of IkeV. A con- 
stant value of Ik eV should suffice for an order of magni- 
tude estimate as iBaeanoff etHI (|2003l ) find that the gas 
temperature varies from approximately 1.9 to 1.3keV from 
rgai = 1-5 to 10" (assuming an optically thin plasma model) . 
iQuataertI ()2004h 's model also predicts that the temperature 
varies from about 2.5 to IkeV from rgai = 0.3 to 10". By 
plugging the val ue of L2-10 within 1.5" (2.4 x lO^^ergs"^) 
as measured by iBaeanoff et al. (2003) into equation (|25|l , 
we find AI ~ 10~^ MQyr~^ . This value is consistent with the 
mass inflow rate at ~ 1" calculated bv lQuataertI l|2004h . 

Our results are not sensitive to the choice of the lower 
limit in the integral across rgai- The lower limit should be 
at most a few 0.01s of pes to at least ~ 10~^pc. The former 
value is the tidal radius for the SMBH at the Galactic centre 
for a IM0 star. Unbound material due to stellar collisions 
or from stellar wind should not exist at smaller radii since 
there are very few stars there to produce it. The value of the 
integral thus ranges from about unity to a few 10s. Since M 
depends upon the square root of this value, the exact value 
of rmin only affects our calculation at the level of a factor 
of a few, and we thus take the square root of the integral to 
be unity. 

Having established that M ~ lO~^M0yr~^ in the vicin- 
ity of 1.5", we now calculate the expected mass loss rates 
due to stellar collisions for different PDMFs. The value of 
M that contributes to the 2-lOkeV ffux is given by: 



M = 



O.OSpc 



dM ^, , ,3 

-73 C,\rTai)d rgai, 

a rgai 



(26) 



where we have shown how to calculate the mass loss rate 
profiles, dM / d^rgai in the previous section. We account for 
the fact that not all of the emission from the unbound gas 
contributes to the 2-lOkeV band with C,{rgai), defined as the 
fraction of fiux from gas at radius rgai with 2keV ^ hv ^ 
lOkeV: 



(.{rgai 



J2kiiV 



dv 



-2keV/kBT(rg^,) 



Jo 



0° Jf 



dv 



-10fceV-/fcBT(rg„,) 



(27) 



Since the gas at each radius is at a slightly different temper- 
ature, and since ( is exponentially sensitive to the temper- 
ature, we must estimate T{rgai)- We do this by setting the 
thermal energy of the unbound material equal to the kinetic 
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Figure 12. The fraction of fiux emitted from unbound material 
at radius Vgai tliat contributes to tlie 2-lOkeV band. 

energy at a radius rgai, and find that 
ksTirgai) ~ mpa\rgai) = 7.8 x lO^'keV (^^^ . (28) 

We plot equation H27|) in Fig. 1121 The value of goes to zero 
at the highest and smallest radii since, for the former, the 
gas is cool and emits most of its radiation redward of 2keV, 
and for the latter, the gas is hot and emits mostly blueward 
of lOkeV. Thus, even though the integral in equation (|26p 
extends to rgai ~ 0, the contribution to M is suppressed 
exponentially at the smallest radii. 

Since, by equations (|20|l and (|2ip . M depends on the pa- 
rameters of the PDMF, we now constrain these parameters 
by limiting the allowed mass loss rate from stellar collisions 
calculated via equation (|26p at 10~^MQyr~^. We consider 
changes in Mmin and a, and keep Mmax set at 125M0 since 
(as seen in Fig.|9]) M is approximately independent of M,nax- 

We sample the Mmin — ce parameter space and use equa- 
tion (|26p to compute the total mass loss rate, the results of 
which are shown in Fig. 1131 The contours represent the cal- 
culated M values, where the solid contours are on a logarith- 
mic scale, and where they are limited from above at a value 
of M = lO~^M0yr~^. The lines are on a linear scale with 
intervals of 1.5 x 10~®MQyr~^, with the thick line denoting 
the 10"^MQyr"^ level. The figure shows that PDMFs with 
flat to canonical-like profiles are allowed. Very top-heavy 
profiles (a < 1.25) are not allowed, as they predict too high 
of a mass loss rate. Mass functions with Mmin ^ TM© are 
also not allowed. These results are consistent with measure- 
ments of the Arches star cluster, a young cluster located 
about 25pc from the Galactic center. Recent measur ements 
l|Kim et al.ll2006l : IStolte et al.ll2005l : iFiger et al.|[l999l ) prob- 
ing stellar masses down to about IMq show that the cluster 
has a flat PDMF, with a in the range of about 1.2 to 1.9 
(depending on the location within the cluster). 

Since M is a much stronger function of a than of Mmin 
it is difficult for us to place tight constraints on the allowed 
range of Mmin- Fig. 1131 shows that we can, however, place 



a constraint on the allowed upper limit of Mmin, since very 
high values of Mmin result mass loss rates > 10~^MQyr~^. 
For a > 1.25, we fit a 3rd degree polynomial (the dashed 
line in Fig. [T3J to the M = lO'^^Moyr"^ contour. This fit 
analytically expresses the upper limit of Mmin as a function 
of a. We provide the coefficients of this fit in the caption of 
Fig.IIl 

The small difference between the solid and dashed lines 
at rgai ~ 0.06pc in Fig. HJd suggests that, even for stellar 
encounters involving small impact parameters, our integra- 
tion does not miss many collisions by ignoring gravitational 
focusing. To estimate the contribution to the total mass loss 
rate in Fig. [13] from gravitational focusing, we take SMtyp, 
the typical amount of mass lost per collision, to be simply a 
function of h. This avoids the multi-dimensional integrations 
involved in equations (|20l and [21) . since for these equations 
Apd is a function of 6, Mpr, Mpd, and Vrei{rgai)- For sim- 
plicity, we choose SMtyp{b) to decrease linearly from 2Mq 
(we assume that both stars are completely destroyed) at 
fe = to at fe = 6o- We find bo by noting from Fig. [Tj 
that for all values of the polytropic index, the amount of 
mass loss for an indirect collision goes to zero at around 
7 — 0.98. By recalling the definition of 7 (equation (O), we 
solve for 60 at 7 = 0.98 by setting Mpr = 1, and taking 
Vret ~ 2a{rgai = 0.06pc). By calculating dF/db (< rgai) (for 
Salpeter values) evaluated at 0.06pc across a range of b, and 
multiplying by 5Mtyp{b), we are able to estimate dM/db. We 
do this for dF/db (< rgai) with and without gravitational fo- 
cusing, and integrate across b. Subtracting the two numbers 
results in our estimate of the contribution to the total mass 
loss rate due to gravitational focusing: 2.3 x lO'^^M©. This 
is about twice the mass loss rate from Fig. [TJ] evaluated at 
Salpeter values. We perform the same calculation across the 
Mmin - ot parameter space, and find that gravitational fo- 
cusing contributes a factor of at most ~ 2.5 to the total mass 
loss rate. 

An underestimate of a factor of 2.5 slightly affects the 
region of parameter space that we are able to rule out, as 
shown by the line contours in Fig. [13] The contours are 
on a linear scale, starting at 4 x lO~^M0yr~^, and ending 
at 2.5 x lO"^M0yr"\ in intervals of 1.5 x lO'^'Meyr'^ 
The 4 X lO"'^M0yr"^ contour (2.5 times less than the lO"'^ 
contour) shows that the region of the parameter space that 
is ruled out is Mmin 1.4M0 and a < 1.4. 

As noted in tj2.3l the amount of mass loss that we cal- 
culate for a single direct collision sometimes over or under- 
pr edicts the amou nt of mass lost (as compared to the work 
of lFreitag fc Bend (2005)) by a factor of a few to at most 
a factor of about 10 (see Figs. [5] and |6]l . We find that our 
calculation of mass loss for direct collisions shows no prefer- 
ence over whether the amount of mass lost is over, or under- 
predicted when considering different combinations of A4\, 
M2, Vrei and b. We therefore suspect that when integrating 
over all of these parameters to obtain the total amount of 
mass lost, our error will roughly cancel. The line contours 
in Fig. [13] serve as a good gauge of how our constraints on 
the allowed region of the parameter change given an uncer- 
tainty of a factor of a few in our calculation of M. As noted 
in the previous paragraph, an underestimation of the total 
mass loss rate results in too high of an allowed upper limit 
in Mmin by a factor of a few, and too low of an allowed lower 
limit in a by about 10%. An overestimation of the total mass 
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Figure 14. a) The IMF power-law slope as a function of the 
PDMF power-law slope for the case of constant star formation, 
b) The same except for exponentially decreasing star formation 
with Texp = 3, 5, 7, 9Gyr (bottom to top line respectively). 



Figure 13. The total mass loss rate contributing the 2-lOkeV 
flux (calculated from equation 126}) as a function of Mmin and 
a. The solid contours are on a logarithmic scale, and are limited 
from above at lO~^M0yr~^. The line contours are on a linear 
scale, and are separated by intervals of 1.5 X 10~®MQyr~^. The 
thick line denotes the 1 X lO~^M0yr~^contour. The dashed line 
is a 3rd order polynomial fit which represents the absolute al- 
lowed upper limit of M^mn as a function of a. The coefficients 
of this polynomial are 21.71, -42.37, 25.33 and -4.27 for ao to 03 
respectively. 



loss rate of a factor of a few results in too high of an allowed 
lower limit in a by a few tens of percents. An overestimation 
will also result in a significantly too low allowed upper limit 
in Mmin since (as seen in the figure) M is a much stronger 
function of a for small a so that the contours are nearly 
vertical. 



6 IMPLICATIONS FOR THE IMF 

We now place constraints on the IMF in the Galactic centre 
with a simple analytical approach that connects the IMF to 
the PDMF, and with the results of the previous section. The 
mass function as a function of time is described by a partial 
differential equation that takes into account the birth rate 
and death rate of stars: 

daM,t) 



dt 



RB{tMM)-^- 



(29) 



'r*(M)' 

where Rsit) is the birth rate density of stars 
(dNB/idtct^rgai)), ${M) is the initial mass function 
normalized such that / <^{M)dM = 1, and Ti,(M) is the 
main sequence lifetime of stars as a function of stellar mass. 
For the initial mass function, we take a power law, 

$ = M-^, (30) 

and for Tt,(M) we use the expression given by 
IMo. van den Bosch, fc White! (|2010l ) 



2.5 X 10^ + 6.7 X 10'^ M^-^ + M*-^ 



3.3 X lO-U'/i-s -H3.5 X 10-iM4 5*'^'''' ^^^^ 
valid for O.OSMq < M < lOOM© and for solar-type metal- 
licity. 

In the following paragraphs, we consider different star 
formation history scenarios. For each scenario, we will need 
to know Rb{tmw), the star formation rate density in the 



Galactic centre at the age of the Milky Way (which we take 
to be 13Gyr). A rough estimate of this value is given by 
the number density of young stars in the Galactic centre 
divided by their age: Rb{tmw) ~ p{r)ri/{{T){M)). Here 
(r) and (Ad) are the average age and average mass of the 
young stars in the Galactic center, which we take to be 
~ lOMyr and ~ lOM© respectively. The parameter 77 is 
the fraction of stars with masses above IOM0, which for 
reasonable mass functions is ~ 0.1%. For self-consistency, 
we use p evaluated at 0.06pc (which from equation (|10p is 
~ lO'^Mopc-^'), since this was the radius with which we 
used to constrain the present-day mass function. These val- 
ues result in Rb{tmw) ~ 10~*pc~'^yr~^. 

For the simple case of a constant star formation rate, 
RB{t) ~ Rb{tmw), and the solution to equation p9)l with 
the boundary condition that ^{M,t = 0) = ${M)ntot{t = 
0), evaluated at the current age of the Milky- way is: 

e(M,t = TMw) = ${M)e-^"'^^^*^^''> 

X {i?s(rMH/)r.(M)e^*"''/"*(*^' 

- RB{TMw)r4M)+ntot{0) | . (32) 

We evaluate the solution at the age of the Milky- way (yield- 
ing the PDMF) because want to compare with our con- 
straints on the PDMF as found in the previous section. 
To solve for ntot{0), we use the known mean density of the 
Galactic centre today at 0.06pc, p{TMw,r = 0.06pc), insert 
equation p2[) into the following expression: 



(33) 



p(rMW, r = 0.06pc) = J £,[M,t = TMw)MdM, 



and solve for ntotiO)- 

We solve for £,[M,tmw) for a range of different IMF 
power-law slopes, 7, and fit a power-law to the solution, 
with a power-law slope a. We plot the IMF power-law slope 
as a function of the calculated PDMF power-law slope for 
constant star formation in Fig. 114b . We have constrained the 
PDMF in the previous section to have a > 1.25, indicated 
in the figure by the vertical line. The figure therefore shows 
that for the case of constant star formation, the IMF power- 
law slope, 7, must be > 0.9. 

For the general case of a star formation rate that varies 
with time, i?s(t) 7^ Rb{tmw), and the solution to equation 
(|29l) with the same boundary condition and evaluated at 
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tmw is: 

£.{M,t = tmw) = $(M)e-"»^"'/"*('^' (34) 

For an exponentially decreasing star formation history, the 
star formation rate is given by: 

RB{t) = i?s(rA/w)e-('-"*^^''' (35) 

Given this star formation history, we solve for £,{M,tmw) 
(by solving for ntot{Q) with equation (I33p ) for Texp =3, 5, 7 
and 9Gyr. We fit power-laws to the resulting PDMFs, and 
show the results in Fig. 114b . The figure shows that smaller 
values of Texp result in larger values of a for any given 7. The 
trend can be explained by the fact that since a smaller value 
of Texp results in a steeper Rb profile, and that all profiles 
must converge to Rb{tmw) at the present-time, Rb pro- 
files with smaller values of Texp have had overall more star 
formation in the past. More overall star formation means 
that the present-day mass function is comprised of fraction- 
ally more lower-mass stars since the IMF favors lower-mass 
stars. The constant build-up of lower-mass stars results in a 
steeper PDMF, so that for any given 7, a should be larger. 
The figure shows that for exponentially decreasing star for- 
mation 7 must be > 0.6, 0.8, 0.9 and 1.0 for Texp = 3, 5, 7, 
and 9Gyr respectively. 

The final case we consider is episodic star formation, 
where each episode lasts for a duration At, where the ending 
and beginning of each episode is separated by a time, T, and 
where the magnitude of each episode is Rb{tmw)- For such 
a star formation history, the solution to equation H34p is: 

^{A4,t = tmw) = $(M)e-^''^»'/^*(*^' 

X \rb{tmw)t4M) [e[(n + l)At + nT]/.,(M) (gg) 

-e"(^'+^'/-('^^'] +n*„,(0) |, (37) 

where n,nax = fioor{(rjv/w — At)/{T + At)}, and where we 
again solve for ntot{0) with equation H33[) . We consider 9 
cases with At and T = lO'', 10^, and 10*^yrs, and show the 
results in Fig. 1151 In each panel the lowest line is At = lO^yrs 
and the highest line is At = lO^yrs. For T = lO^yrs, 7 > 
0.8 and 0.5 for At = 10^ and lO'^yrs respectively, while the 
At = lO^yrs case results in constraints on 7 that are too low 
to be realistic. For T = lO'^yrs, 7 > 0.5 and 0.4 for At = 10^ 
and lO^yrs respectively, while again, the At = lO^yrs case 
results in unrealistic constraints. Finally, for the T = lO^yrs, 
7 > 0.5 for At = 10", while the At = lO'^ and lO^yrs case 
result in unrealistic constraints. We test whether when the 
last star formation episode occurs relative to the present day 
affects our solution of £,{M, tmw) by varying the start time 
of the star formation episodes. By varying the start time and 
testing all the combinations of At and T that we consider, 
we find that the lines in Fig. [15] vary by at most about 5%, 
so that the main trends in the figure are unaffected. 
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Figure 15. The IMF power-law slope as a function of the PDMF 
power-law slope for the case of episodic star formation. In each 
panel the lowest line is At = lO^yrs and the highest line is At = 
10*^ yrs. 



7 CONTRIBUTION FROM RED GIANTS 

Spectroscopic observations have revealed that the central 
parsec of the Galaxy harbors a significant populat ion of gi- 
ant stars (|Paumard et al.ll200^ : iManess et al.l^2007^ . Due to 
their large radii (and hence large cross sections), it is possi- 
ble that they could play an important part in the mass loss 
rate due to collisions in the Galactic centre. 

In assessing their contribution to the mass loss rate, 
care must be taken when deriving the collision rat es, because 
their radii, rue, are strong functions of time, t. iDale et al.l 
(2009) have already calculated the probability, P{rgai), for 
a red giant (RG) in the Galactic centre to undergo collisions 
with main sequence impactors. They have taken into account 
that Taait) by integrating the collision probability over the 
time that the star resides on the RG branch. We use their 
results to estimate the mass loss rate due to RG - MS star 
collisions. 

To find the number density of RGs in the Galactic cen- 
tre, we weight the total stellar density by the fraction of 
time the star spends on the RG branch: 

nRcirgai) ^ n-i,{rgai)^-^- (38) 
r* 

This approximation should be valid given a star formation 
history that is approximately constant when averaged over 
time periods of order thq. The number of collisions per unit 
time suffered by any one red giant, P{rgai), should of order 
the collision rate averaged over the lifetime of the RG, and 
is given by 

PM-(PM) (39) 

\ It TRG 

If we define 5M to be the typical amount of mass lost in the 
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collision, then the mass loss rate is 



dM 



dlnr 



= 47rr 



dM 



gal 



gal 



d^r, 



gal 



gal) 



gal, 



TRG 



SM. 



(40) 

To calculate an upper limit for the contribution of RG 
- MS star collisions to the mass loss rate, we assume all 
RG and MS stars have masses of IM©, and that the en- 
tire RG is destroyed in the collision. Collisions involving 
IM0 RGs yield an upper limit, because there is not an 
appreciable amount of RGs with masses less than ~ IMq 
due to their MS lifetimes being greater than the age of the 
Galaxy. For RGs with masses greater than IM© , the amount 
they contribute to the mass loss rate is a competition be- 
tween their lifetimes and radii. Red giant lifetimes decrease 
with mass (thereby decreasing the time they have to col- 
lide) and their radii increase with mass (ther e by in creasing 
the cross section). In their Fig. 3. iDale et al.l (|2009l ) clearly 
show that the number of collisions decreases with increas- 
ing RG mass, indicating that the brevity of their lifetime 
wins over their large sizes. One solar mass MS impactors 
should yield approximately an upper limit to the mass loss 
rate, since ~ IM© MS stars are the most common for the 
PDMFs under consideration. 

Since we assume that the entire RG is destroyed in the 
collision SM — IMq. For the case that all impactors are 
IM0 MS stars, we calculate URcirgai) from equation (|38p by 
noting that ni,{rgai) ~ p-kifgai) / (IMq). For self-consistency, 
we must truncate P{rgai) at 1 for all Pifgai) > 1 since we are 
considering the case where one collision destroys the entire 
star. We plot equation H40|) for this calculation in Fig. 1161 
The discontinuity is due to our truncating P(rgai) at 1. The 
figure shows that the mass loss rate for RG-MS star collisions 
never exceeds lO~^M0yr~^, well below typical dM /dlnrgai 
for values for MS - MS collisions (see Figs . l9l and [TOll . More- 
over, in their hydrodynamic simulations, iDale et all (|2009l ) 
note that in a typical RG - MS star collision, at most ~ 10% 
of the RG envelope is lost to the RG. We therefore conclude 
that the contribution of RGs to the total mass loss rate in 
the central parsec of the Galaxy is negligible. 

The figure shows that by Vgai ~ 0.06pc, the mass loss 
rate for RG-MS star collisions is at most about lO~^M0yr~^. 
It is thus possible that for MS-MS collisions, values of 
Mmin and a that results in total mass loss rates just be- 
low lO~^M0yr~^ could be pushed past this threshold with 
the addition of mass loss due to RG collisions. However, we 
believe that this is unlikely for two reasons. The inclusion of 
the factor, C^, when calculating the total mass loss rate (see 
equation (|26p ') will reduce the mass loss by at least a factor 
of 0.6 (see Fi g. 1121). Also, as n oted by the hydrodynamic 
simulations of lDale et al. M2OO9I ). for a typical RG - MS star 
collision, at most ~ 10% of the RG envelope is lost to the 
RG. This will reduce dM /dlnrgai for RG-MS collisions by 
another factor of 10. 



8 CONCLUSIONS 

We have have derived novel, analytical methods for calcu- 
lating the amount of mass loss from indirect and direct stel- 
lar collisions in the Galactic centre. Our methods compares 
very well to hydrodynamic simulations, and do not require 
costly amounts of computation time. We have also computed 




Figure 16. An upper limit to the mass loss rate due to collisions 
between RG and MS stars. The arrow indicates the range in the 
x-ray observations (rg^j < 1.5") with which we use to constrain 
the PDMF (see 



the total mass loss rate in the Galactic centre due to stel- 
lar collisions. Mass loss from direct collisions dominates at 
Galactic radii below ~ O.lpc, and thereafter indirect colli- 
sions dominate the total mass loss rate. Since the amount 
of stellar material lost in the collision depends upon the 
masses of the colliding stars, the total mass loss rate de- 
pends upon the PDMF. We find that the calculated mass 
loss rate is sensitive to the PDMF used, and can therefore 
be used to constrain the PDMF in the Galactic centre. As 
summarized by Fig. 1131 our calculations rule out a < 1.25 
and Mmin TM© in the Mmin — ct parameter space. Finally, 
we have used our constraints on the PDMF in the Galactic 
centre to constrain the IMF to have a power-law slope > 
0.4 to 0.9 depending on the star formation history of the 
Galactic centre. 
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